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CN ■ ABSTRACT 



Context. LHS 1070 is a nearby multiple system of low mass stars. It is an important source of information for probing the low mass 
end of the main sequence, down to the hydrogen-burning limit. The primary of the system is a mid-M dwarf and two components 
are late-M to early L dwarfs, at the star-brown dwarf transition. Hence LHS 1070 is a valuable object to understand the onset of dust 



formation in cool stellar atmospheres. 

Aims. This work aims at determining the fundamental stellar parameters of LHS 1070 and to test recent model atmospheres: BT-Dusty, 
BT-Settl, DRIFT, and MARCS models. 
Methods. Unlike in previous studies, we have performed a ^-minimization comparing well calibrated optical and infrared spectra 
with recent cool star synthetic spectra leading to the determination of the physical stellar parameters T c g , radius, and log g for each of 
C/j , the three components of LHS 1070. 

Results. With exception of the MARCS models which do not include dust formation, the models are able to reproduce the observations 
and describe the main features of the visible to IR spectra. This is consistent with the fact that dust formation prevails in the B and 
C component atmospheres. The parameters obtained with the DRIFT models confirm the values determined in earlier studies. But 
important differences between models are observed, where the MARCS model is too bright in the H and K bands, and the BT-Settl 
and BT-Dusty models systematically yield up to 100 K higher r c ff in the case of the B and C components. This confirms a trend 
for models without, or with less efficient cloud formation, to predict higher r cff than models richer in dust (DRIFT). Even models 
including cloud physics however still produce slightly too bright J band flux, showing as too blue J - K colors. The onset of dust 
formation remains therefore a particularly challenging regime to understand. 
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1. Introduction tant role in understanding the formation of brown dwarfs and 

planets. 

The lower end of the Hertzsprung-Russel diagram has much 

importance as the vast majority of stars in the Galaxy are low Despite their large number in the Galaxy, little is known 

mass stars. In the Galaxy, 70% of the stars are M dwarfs, about low mass stars because of the difficulty (i) to get a ho- 

c3 \ They contribute o v er 40% of the t o tal stellar m ass content mogeneous sample with respect to the age and metallicity due 

dGould et all Il996t iMera et all Il996t iHenrvl Q998). These M to their intrinsic faintness, and (ii) to disentangle the parameter 

dwarfs have a mass that ranges from 0.6M o to the hydrogen space (T eS , logg, and metallicity). Indeed a number of studies 

burning lim it of about 0.075 to 0.0 85Mq depending on the have shown that a change in temperature or gravity can com- 

metallicity (IChabrier & Baraffd |2000|) . These stars are found pensate for a change in metallicity to some degree. An addi- 

in any p opulation, from young metal r ich M-dwarfs in open tional difficulty is the complexity of their atmospheres: convec- 

clusters dReidl 11993b lLeggettetal1ll994l) to the several billion tion in optically thin regimes, molecules, and dust cloud forma- 

years old metal poor dwarfs in the galactic halo (Green and tion for the later types. Water vapor and CO bands dominate 

Morgan 1994) and in the globular clusters dCool et al.lll996t the Rayleigh-Jeans branch of the spectral energy distribution 

iRenzmiet ail 19961) . Such low mass stars are an important probe at infrared wavelengths (> 1 .3 fan), while TiO, VO, and metal 

for our Galaxy as they carry fundamental information regard- hydrides govern the corresponding visual (> 4000A) to near- 

ing the stellar physics, galactic structure and formation, and its infrared (< 1.3 //m) spectral energy distribution. Convection 

dynamics. In addition, the existence of brown dwar fs or planets reaches out to the optically thin (as far as to an optical depth t ~ 

being discovered and con firmed around M-dwarfs ( Butle r et al.l 10 3 ) portion of the atm ospheres, flattenin g the temperature gra- 

20041 iBonfils et al.l201 ll and references therein) plays an impor- dient of the atmosphere (Allar d et al.ll997h . lLudwig et alJ (f2002. 
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2006) have determined the mixing length based on a compari- 
son of the mixing length theory used in ID static model atmo- 
spheres and Radiation Hydrodynamic simulations. In M dwarfs 
later than M6 the outermost temperatures fall below the con- 
densation temperatures of silic ate grains, which l ea ds to the for- 



mation of dust clouds (see e.g. Tsuii et alj[l 996a,b: Allar d et al 



1997; iRuiz et al J [19971 [Allard & Hauschildtl Tl998l: lAllard et al 



1998). These processes complicate the understanding of these 
cool atmospheres. 

One approach to study the physics at the low end of the main 
sequence is to compare observed spectra with synthetic spectra 
from various authors and modelling techniques. The determina- 
tion of the physical parameters (effective temperature, gravity, 
metallicity, radius) is obtained by spectral synthesis, i.e. x 2 min- 
imization. 

LHS 10 70 is a low mass m ultiple system of cool dwarfs dis- 
covered bv lLeinert et al.l (1 19941) . with visual magnitu de 15. It is 
locate d at a distance of 7.72+0.15 pc from the Sun ( Costa et alj 
120051) and is considered as a memb er of the disk population 
with a probable age of around 1 Gyr dReiners et al. 2007b). The 
spectral types for the A, B, and C components were found to 
be M 5.5-M6, M8.5, and M9-M9.5, respectively dLeinert et alJ 
2000). A fourth component was suspected very close to the pri- 
mary bv lHenrv et al.l (1 19991) from HST Fine Guidance observa- 
tions, but this detection is no longer considered to be real (T. 
Henry, private communi cation). The l atest orbit determination 
has been performed by iKohler et al.l d2012l) . with semi major 
axes of 0.458" for the close pair BC and 1.112" for the wide 
orbit of BC around component A. E.g., on December 12, 2003, 
component B was separated from A by 1.77" at position 13°, and 
co mponent C from B by 0.41" at 178°. 

iLeinert et al.l d!998l [2000) have derived effective tempera- 
tures of 2950 K, 2400 K and 2300 K for the components based on 
spectral analysis, and found that B and C showed clear signatures 
of dust in their spectra. They presented photometric mass esti- 
mates ranging from 0.109 to O.O79M for the three stars based 
on theoretical isochrones, thus reaching right down to the mini- 
mum hydrogen burning mass. This mass range makes LHS 1070 
a valuable system for understanding the formation of dust in cool 
atmospheres and the processes that occur at the star/brown dwarf 
transition. LHS 1070 is therefore a testbed to validate and define 
further developments of both atmospheric and interior models 
at the lower end of the main sequence. We assume the same age 
and composition for the three components of this system for sim- 
plicity. 

In this paper, we present the spectral synthesis of the 
components A, B, and C. We determine their physical pa- 
rameters by comparing the well-calibrated HST spectra in 
the optical (from FOS) and in near and mid-IR (from 
HST/NICMOS and ISOPHOT-S) with synthetic spectra com- 
puted fro m recent stellar atmosp here model s: BT-Dusty and 
BT-S ettl dAllard & Homeierll2012l). MARCS dGustafsson et alJ 
120081) . and DRIFT dWitte et alj|2009l) . Observations and data re- 
duction are described in Sec. [2] Sec. [3] presents the atmosphere 
models used in the analysis. In Sec.[4]we give the determination 
of the stellar parameters and show the comparison between ob- 
served and modeled spectra. Discussion and conclusion follow 
in Sect. [5] 

2. Observations and data reduction 

Regarding th e optical photometry and spectroscopy, the reader 
is referred to ILeinert et alJ |[T000). The new input concerns the 
infrared data. 



2.1. Photometry 

The J, H, K, and L phot ometry pres e nted i n Table 1 for 
LHS 1070 ABC refers to iLeggett et al.l (fl998h . We obtained 
M-band photometry with the UIST instrument on UKIRT on 
November 9, 2002 , as well as N-ban d photometry with the 
MAX camera ( Robb erto & Herbsdfl998l) on the same telescope 
on August 27, 1996. The brightnesses of the individual com- 
ponents wer e then derived from the brightness rati os obtained 
with NACO dRousset et all200llLenzen et al.ll200l for J, H, K 
on December 12, 2003, for L and M on December 6, 2001. The 
MAX observations in N-band yielded separate brightnesses only 
for component A and the sum B+C of the other two components. 

The M band data were obtained in the Mauna Kea 
Observatory Near-Infrared System (/l cen trai = 4.7/an, 50%- 
width= 0.23/mi) in service mode. Aperture photometry was ap- 
plied to the pipeline-reduced data. The N band photometry used 
a standard filter (/l ce ntrai = 10.47/zm, 5%-width = 4.65/vm). After 
the standard processing steps (bad pixel correction, combina- 
tion of individual chop cycles after shift-and-add) we performed 
aperture photometry. The absolute calibration at 10.4 fi relied on 
the HR 6464 spectrophotometric standard mo del, produced by P. 
Hamm ersley using the procedure described in lHammerslev et aLl 
(Il998l) . and made available on the ISO web page for ISO stan- 
dards. HR 6464 has the same spectral type MOM as the stan- 
dard HR 400 actually observed, and the flux ratio of the two was 
determined from their fluxes in the IRAS 12 fim band. Strictly 
speaking, the result is a narrow-band (0.25 fim) brightness at 
10.4 fim under the assumption of an MOM type spectrum. The 
spectral slopes in this wavelength region are smooth and all are 
representing Rayleigh- Jeans tail emission. Uncertainties result- 
ing from the difference in spectral type between MOM and our 
object therefore are not important. To be conservative, errors of 
+ 5 mJy, + 4 mJy, and + 3mJy are taken for the fluxes of the 
combined system, of component A, and the sum of components 
B and C, respectively. 

Our NACO observations were used in determining the rela- 
tive brightnesses of the three components of LHS 1070, because 
of the superior spatial resolution of this instrument. After stan- 
dard reduction (flat fielding, bad pixel correction), aperture pho- 
tometry was applied. The NACO bandpasses closely match the 
MKO-NIR system for the J, H, K, L bands. The main difference 
is in the N band (4.8 fim for NACO versus 4.7 fim at Mauna Kea, 
20% width of 0.64 /mi for NACO versus 50% width of 0.23 fim 
at Mauna Kea). In this wavelength range, the slopes of the com- 
ponent spectra are very similar. Therefore the NACO-measured 
brightness ratios were applied to extract the component bright- 
nesses from the MKO-NIR based total brightnesses. 

LHS 1070 was observed with ISOCAM (ICesarskv et al.l 
[l996h on the ISO satellite in the LW2 (6.7 fim) and LW3 (14.3 
fim) filters on November 28, 1996 (PI H. Zinnecker). Near- 
infrared narrow and medium-band photometry from 0.90 fim to 
2.15 fim was obtained with the NIC1 and NIC2 cameras of the 
HST NICMOS instrument on Jan. 2, 1998. For the analysis we 
used pipeline reduced images. On each frame the system was 
clearly resolved, and it was possible to obtain separate photo- 
metric measurements for the three stars. Before performing aper- 
ture photometry for a given star, we removed the images of the 
other two objects by subtracting a scaled PSFPSF template, con- 
structed from component A and shifted to the known coordinates 
of the stars. An aperture correction, taken from the NICMOS 
Data Handbook, was applied, and photometric calibration was 
performed via multiplying by a conversion factor between the 
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Table 1. Photometric data. Fluxes are in logl0(ergs crrT 2 s A~'). 



Wavelength 


Filter 1 


Component A 


Component B 


Component C 


Components A+B+C 


//m 








logl0(FJ 




logl0(F,i) 


original notation 


0.900 


NIC1 F090M 


-13.347 ± 


0.010 


-14.122 + 0.016 


-14.289 ±0.017 




156 ± 2 mJy 2 


0.953 


NIC1 F095N 


-13.213 ± 


0.031 


-13.870 ±0.03 


-14.056 ±0.036 




253 ± 8 mJy 2 


0.970 


NIC1 F097N 


-13.251 ± 


0.031 


-13.943 + 0.027 


-14.092 ±0.028 




237 ± 7 mJy 2 


1.083 


NIC1 F108N 


-13.229 ± 


0.024 


-13.831 ±0.017 


-13.953 ±0.017 




332 ± 7 mJy 2 


1.100 


NIC1 F110M 


-13.292 ± 


0.006 


-13.887 ±0.007 


-14.037 ±0.009 




295 ± 2 mJy 2 


1.130 


NIC3 F113N 


-13.267 ± 


0.022 


-13.809 ±0.073 


-13.965 ±0.097 




343 ± 15 mJy 2 
355 ± 7 mJy 2 


1.130 


NIC1 F113N 


-13.249 ± 


0.020 


-13.796 ± 0.019 


-13.962 ±0.022 




1.25 


J 


-13.350 ± 


0.017 


-13.904 ±0.028 


-14.036 ±0.028 


-13.189 + 0.012 


9.14 ±0.03 mag 


1.450 


NIC1 F145M-13.494 ± 0.005 


-14.077 ± 


0.005 


-14.218 ± 0.005 




326 ± 2 mJy 2 




1.640 


NIC1 F165M 


-13.517 ± 


0.004 


-14.026 ± 0.004 


-14.163 ±0.004 




419 ± 2 mJy 2 


1.65 


H 


-13.535 ± 


-0.018 


-14.061 ± 0.028 


-14.195 ±0.031 


-13.367 ±0.013 


8.51 ±0.03 mag 


1.660 


NIC3 F166N 


-13.492 ± 


0.013 


-13.930 ±0.084 


-14.130 ±0.133 




472 ± 29 mJy 2 


1.800 


NIC2F180M 


-13.677 ± 


0.006 


-14.197 ±0.008 


-14.333 ±0.009 




346 ± 2 mJy 2 


1.900 


NIC1 F190N 


-13.751 ± 


0.014 


-14.276 ±0.007 


-14.415 ±0.008 




324 ± 4 mJy 2 
314±2mJy 2 
378 ± 5 mJy 2 


2.040 


NIC2 F204M 


-13.833 ± 


0.006 


-14.343 ±0.004 


-14.474 ±0.004 




2.150 


NIC2F215N 


-13.825 ± 


0.013 


-14.273 ±0.005 


-14.400 ±0.006 




2.2 


K 


-13.860 ± 


0.018 


-14.339 ± 0.018 


-14.468 ±0.024 


-13.674 ±0.014 


8.14 ±0.03 mag 


2.300 


NIC2 F222M 


-13.844 ± 


0.004 


-14.301 ±0.003 


-14.429 ±0.003 




378 ± 2 mJy 2 


2.375 


NIC2 F237M 


-13.973 ± 


0.005 


-14.447 ±0.003 


-14.579 ±0.003 




317±2mJy 2 


3.8 


L' 


-14.604 ± 


0.013 


-14.990 ±0.013 


-15.097 ±0.014 


-14.377 ±0.012 


7.63 ± 0.06 mag 


4.78 


M 


-15.008 ± 


0.019 


-15.449 ±0.028 


-15.536 ±0.041 


-14.806 ±0.017 


7.72 ± 0.04 mag 


6.7 


LW2 










-15.349 ±0.019 


70 ± 3 mJy 


10.4 


N 


-16.190 ± 


0.053 


-16.378 


± 0.079 


-15.973 ±0.035 


38.4 ± 3 mJy 


14.3 


LW3 










-16.555 ±0.043 


21+2 mJy 



1 M: medium-band filter (AA = 0. 1 - 0.2/mi), N: narrow-band filter (AA = 0.02 - 0.04/mi) 

2 these values are the sum of the values measured directly for the individual components 



counts and F v as stated in the Handbook. The results are given 
in Table 1 . 

2.2. Spectroscopic observations 

Near-infrared spectra of the individual components of LHS 1070 
were taken with grisms G096, G141 and G206 and the NIC3 
camera of the HST NICMOS instrument on January 2, 1998 
for the J, H, and K bands. In a first step we removed from the 
four spectra we had for each grism the short wavelength and 
long wavelength ends. For grism G096 this left the range of 
0.80-1.085 /mi or 0.8-1.15 /mi, depending on the local noise 
level. Similarly, G141 covered the range of 1.10-1.60 /mi or 
1.10-1.85 /mi. G206 ranged from 1.65 /mi to 2.45 /mi. For 
each wavelength pixel, all spectra covering that wavelength 
were averaged using a weighted average, where the weights 
came from the formal uncertainties of the spectra (typically 
from 2% to 5%). The reduction was based on the NICMOSlook 
IDL-based data reduction package. For the extraction of the 
spectra, the 'no weighting' option was used. The result for the 
A component was slightly scaled (by a few percent) to fit the 
NICMOS photometry. The extraction procedure also gave the 
combined spectrum of components B and C (separated only 
0.4" on the detector), which was decomposed into the spectra 
of the individual components on basis of the brightness ratio 
as a function of wavelength as determined from the narrow- 
and medium band photometry. The spectral resolution of the 
NICMOS spectra is R = A/AA = 200. 

Higher resolution spectra with R = 400, 1500, anc/1400 in 
the J, H, and K ba nds, respectively, w e re taken with the V LT 
NACO instrument dRousset et al J 120031: ILenzen et al.ll2003l) on 
telescope UT4 of the VLT for the individual components of 
LHS 1070. The camera/grism/filter settings were S27/Grism4/J, 



S54/Grism3/H, and S27/Grism3/K, respectively. The slit width 
was 86 milli-arcsec in all three bands. After subtracting the sky 
from the raw data, a flat normalized along the dispersion direc- 
tion was applied. The spectra were then traced and extracted. 
Finally, the individual exposures were averaged. This procedure 
was also applied to the telluric standards. All telluric standards 
are of spectral type G2V, which allows a proper modelling 
of their intrinsic spectrum with that of the Sun. The Solar 
sp ectrum was constructe d using the data and scripts provided 
bv lMaiolino et alJ (U996). Wavelength calibration was obtained 
by arc lamp exposures. The additional effort to guarantee the 
absolute flux calibration was not taken, so only the relative 
shape of the spectrum was determined in each band. For each 
component, we used one scaling factor for each band to bring 
the spectra to the absolute level as resulting from the NICMOS 
observations. The errors, taken as the formal uncertainties of the 
averaging process, typically range from 5-10%. The number of 
pixels per resolution element in the J, H, and K band is about 
3, 1.5, and 3. Because of the missing independent absolute 
calibration, the NACO spectra are used primarily for the study of 
spectral features. However, comparison with the NICMOS spec- 
tra, like in FigureQ] gives confidence also in their spectral shape. 

Spectra for the system from 2.5 /mi to 1 1 .7 //m were obtained 
with the ISOPHOT-S spectrometer on board the ISO satellite 
on November 23, 1996 (PI: T. Tsuji) with an exposure time 
of 4096 s. ISOPH OT-S was the spectroscopic sub-instrument 
of ISOPHOT dLemke et alJ [l996). It had an entrance aperture 
of 24"x24", covering simultaneously the 2.5 to 4.9 and 5.9 
to 11.7 /mi ranges. The spectral resolution of both channels 
ranged between 65 and 130. We processed the ISOPHOT-S 
observa tions using the Pho t Interactive Analysis (PIA) 
V10.S dGabriel et alj|1997l) following the standard data reduc- 
tion scheme. The measurements were further reduced following 
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LHS1070 A NICMOS/NaCo comparison 



LHS1070 B NICMOS/NaCo comparison 



LHS1070 C NICMOS/NaCo comparison 




1.5 2.G 

wavelength (micron) 




1.5 2.0 
wavelength (micron) 




1.5 2.0 
wavelength (micron) 



Fig. 1. Comparison between NICMOS and NACO spectra of the three components. 



our self-developed processing scheme (Kospal et al. submit- 
ted), correcting for the slight off-centre positioning of the source. 



2.3. Spectroscopic features 

Figure [2] shows the optical spectra of all the three components, 
with expected atomic and molecular features in the optical range 
between 5000 to 8500A. The most important ones are molecular 
bands of TiO, CaH, and VO, and atomic lines like Cal, Nal, and 
KI. The TiO bands get weaker towards lower temperature, from 
component A to C, due to condensation into dust species. The 
hydride CaH at 6380 A and 6880 A decreases in strength with 
decreasing temperature. The KI doublet at 7665 and 7699 A 
is very strong and is useful for gravity determination. The Nal 
doublet at 8183 A and 8195A is strong in all the components, 
whereas Cal at 6103 A is weak at all temperatures. 

Figur^!] shows the near-IR spectra obtained with NICMOS 
on HST for the three components as well as the thermal in- 
frared spectrum of the unresolved system A+B+C taken with 
ISOPHOT on ISO. Photometric results are superimposed on 
these spectra. 

Figure [4] shows the NACO spectra obtained at the VLT of 
the three individual components in the / (upper panel), H (mid- 
dle panel), and K (lower panel) band for components A, B, and 
C (from left to right). The main expected features are also indi- 
cated. The Paschen Beta and Cal lines can be seen in the /-band 
spectra with "equivalent widths" of 2-2.5 A and 1 .5-2 A, respec- 
tively. The //-band spectra contain many relatively weak absorp- 
tion features, which defy definite identifications, with possible 
exception of Mg (1.711 /mi), OH (1.689 /mi), and Al (1.675 
/mi). H2O bands define the shape of the / and H band peaks. 
Water absorption is most obvious in the /-band at 1.33 /mi and 
strengthens through the later types: the flux ratio between the 
peak and the minimum of the absorption band increases from 
1.09+0.01 for component A to 1.30+0.01 and 1.39+0.01 for the 
cooler components. The Zf-band spectra of the three components 
show strong CO bands and more or less pronounced atomic fea- 
tures. The Nal lines weaken from the hotter component A to the 
cooler components B and C as dust forms. 



3. Model atmospheres 

The atmospheres of low-mass stars are mostly composed of 
molecular hydrogen and CO, and their spectra characterized by 
strong TiO bands in the optical, water vapor bands in the in- 
frared, chromospheric activity, flares, magnetic spots, and planet 



like properties (lAllard et alj 119971) . The discovery of the for- 
mation of dust clouds in M and L dwarfs makes it even more 
challenging to understand their spectral properties. An effec- 
tive temperature of 2600 K and below is sufficiently low to 
give rise to enough silicate dust formation in the photospheric 
layers to affect the spectral properties of late type M dwarfs. 
These grains produce a "veiling" in the optical by dust scatter- 
ing and an important greenhouse effect (redistribution of flux 
to the infrared ) which strongly i nfluence the infrare d spec- 
tral properties. iTsuii et alj d 19991) : lAllard et all d2001l) treated 
the dust formation in the low-mass stars in pure Chemical 
Equilibrium (hereafter CE). CE yields the formation of con- 
densates from zirconium dioxide and silicates (Mg2Si04) to re- 
fractory ceramics (CaTi0 3 , A1 2 3 ), salts (CsCl, RbCl, NaCl), 
and ices (H2O, NH3, NH4SH4) depending on the tempera- 
ture of the atmosphere from M through T a n d Y spectral 
types dAllard et alJ200ltlLodders & Fee lev 20061. lHelling et al.1 
( 2008ailbh: lAllard & Homeierl d2012l) explored the properties of 
the formation of the dust clouds in low-mass stars by determin- 
ing the radial distribution and average size of the grains. 

We compared the flux- calibrated spectroscopic data of 
LHS 1070 with synthetic spectra computed from three recent at- 
mosphere models for cool stars. These models are described be- 
low. In order to match the observed spectra, the synthetic spectra 
have been scaled by the dilut ion factor (R/d ) 2 where d is the dis- 
tance of the system from Costa et al. (2005) and the stellar radius 
R is a free parameter ranging from 0.096 R e to 0. 142 R Q at a step 
of 0.002. The radius resulting from the fits can then be compared 
to predictions from stellar evolutionary models. 

3. 1 . BT-Dusty and BT-Settl 

Recently, the preliminary res ults of the BT- C ond/ Dusty and BT- 
Settl models based on the lAsplund et al.l d2009t) solar abun- 
dances have been published in a review by Allard & Homeierl 
d2012h . and those grids distributed via the PHOENIX web simula- 
tor ( [http://phoenix.ens-lyon.fr/Grids/BT-Settl/AG SS2009/). For 
this paper, we have developed a revision to these BT-Dusty and 
BT-Settl models based on the ICaffau et al.l d201 lh solar abun- 
dances using slightly revised atomic and molecular opacities as 
well as cloud physics. Their detailed publication is in prepara- 
tion. We summarize below the differences betwee n the current 
versio n of the models and t he published versions of lAllard et al.l 
d200ll) : lAllard & Homeierl d2012l) . 

The BT appellation stands fo r the Barber and Te nnyson so- 
called BT2 water vapor line list (Barb er et alJ l2006). The mod- 
els are provided in several versions addressing different limit- 
ing treatment of the c loud physics . The D usty and Cond appel- 
lation, as defined by lAllard et all d200ll) . refer to a CE treat- 
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Fig. 2. Optical Spectra of the components of LHS 1070 obtained with the Faint Object Spectrograph (FOS) on HST. The atomic and 
molecular features visible in all three components are shown in the upper panel. 



ment of dust formation where dust opacities (assuming spher- 
ical grains with interstellar grain size distribution) are ignored 
in the Cond models. The BT-Dusty and BT-Cond models re- 
produce the color properties of the 2001 AMES Dusty and 
Cond models in the brown dwarf regime, with minor differ- 
ences relating to revised opacities for H2O, TiO dPlez| [l998). 
VO and most o f the hydrides, d etailed damping profiles for 
the alkali lines (lAllard et al.1 120071) . the damping constants of 
molecular lines dHomeier 20Q5| ), new line lists for C H4 (STDS, 
lHomeiedl2005h CO? (CDSD. iTashkun et alj Hjog), and CIA 
dBorvsow et al.ll200Tt I Abel et al.ll2.01 ll) to mention the most im- 
portant chan£es] ; _Were^ervethe NextGen appellation, as chosen 
for the HauschildtetalJ d 1999b models, for corresponding pure 
gas phase models. 



The BT-Settl models described bv lAllard & Homeierld2012l) , 
on the other hand, include gravitational settling which is ig- 
nored in Dusty models, and which involves a cloud model. The 
growth of the grains is governed by the supersaturation ratio of 
the gas and is triggered by the actual collision rates between 
grains and gas molecules and, therefore, depends on the grain 
surface. A brea kthrough was therefore achieved, compared to 
earlier versions d Allard et all I2003L l2007t iReiners etaT]l2007at 
iHelling et alJl2008ch . by calculating locally the supersaturation 
instead of assu ming the fixed conservative value proposed by 
iRossowl dl978l). The clo ud model used in the BT-Settl models is 
based on lRossowl £1978) and accounts for nucleation, condensa- 
tion, supersaturation, gravitational settling or sedimentation, and 
advective mixing. The latter includes convective mixing based 
on the mixing length theory in the convect i on zon e, an exponen- 
tial overshoot according to Ludwig et al. (2006), and the grav- 
ity waves according to lFrevtag et al.l d2010h . The latter dominate 
the mixing in the uppermost atmosphere layers, i. e. the regions 
where clouds first begin to form in late M dwarfs. The grains are 
still considered spherical but the grain sizes (a unique or mean 
value per atmospheric layer) and grain number densities are de- 
termined by the comparison of the different timescales, and thus 
vary with depth to reach grain sizes of a few times the interstel- 



lar values at the cloud base. The BT-Settl models do not enforce 
grains to be in equilibrium with the gas phase as is the case of 
the Dusty and Cond models. However, the gas phase opacities 
reflects the depletion of elements from the gas phase caused by 
grain growth. The numerical solving on-the-fly of the gas phase 
CE in the BT-Settl models (unlike what is often done by other 
authors) costs computing time but allows to account for the cool- 
ing history of the atmospheric layers. 

For the current BT-Settl model version, we have addition- 
ally explored the effects of nucleation in our cloud model by 
assuming a constant nucleation rate of 1 cirT 3 . This acts as a 
limiting factor to the creation, sedimentation and depletion of 
grains and thus leads to the persistence of grains higher up in 
the atmospheric structure than assuming unlimited grain pro- 
duction. We are using in general atomic damping constants ac- 
cording to lUnsoldl dl968[) with a correction fact or to the van der 
Waals widths of 2.5 ( Valenti & Piskunovll 19961) . van der Waals 
broadening of molec ular lines with generic widths according to 
iHomeier et alJd2003l) . More accurate broade ning data for neutral 
hydrogen collisions bv lBarklem et al.l (|2000) have been included 
for several important atomic transitions such as the Alkali, the 
Ca 1 and Ca n resonance lines. 

Grains are assumed here to be spherical and non-porous, 
and their Rayleigh to Mie scattering and absorptive properties 
are considered for 55 types of mater ial including the sp ecies 
mentioned above an d presented in All ard et ail d200ll) and 
iFerguson et al.l d2005l) . plus MnSi03 and Na20sSi2- The grain 
opacities are computed in each layer for pure spherical grains 
using the radius determined by the cloud model (assuming it 
represents an average radius of the grains in each layer) and the 
material density as measured in the laboratory. The opacity con- 
tributions of the various grain species present in each layer are 
finally summed. 

The model atmospheres and synthetic spec tra are com- 
puted with the PHOENIX radiative transfer c ode dAllardl [l990: 
Allard & Hauschildtll 19951 lAllard et al.ll200lb using hydrostatic 
equilibrium, convection based on the mixing length theory and a 
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Fig. 3. Near-IR spectra (solid line) and photometry (solid circles) obtained with NICMOS on HST for components A (upper left), 
B (upper right), and C (lower left). The ISOPHOT thermal infrared spectra of the unresolved system is shown on the lower right 
panel, again with photometric measurements overplotted. 



mixing length which varies from 2.2 to 1.6 from brown dwarfs 
to the Sun according to results of radiation hydrodynamical sim- 
ulations dLudwig et alJ ll999. 2002L 120061) . spherically symmet- 
ric radiative transfer using radii provided by published evolu- 
tion models, micro-turb ulence velocities fr om radiation hydro- 
dynamical simulations ( Frevtag et al.1 12010). and the latest solar 
abundances bv lCaffau et al.l d2011l) . 

The synthetic spectra are provided over the entire spectral 
range of interest at a spectral resolution of 0.05 Ain the optical 
and 0. 1 Ain the near-IR. For this paper, we explored these model 
grids with parameters described as follows: 

- r eff from 2000 K to 3 100 K with 100 K step, as expected for 
mid-M to L dwarfs, 

- logg =4.5,5.0, and 5.5 dex, 

- [M/H] = -1.5,-1.0, -0.5,0.0, +0.3,and+0.5dex, 

3.2. MARCS 



The MARCS code (Gustafsson et al. 2008) assumes hydrostatic 
equilibrium, Local Thermodynamic Equilibrium (LTE), chemi- 
cal equilibrium, homogeneous plane-parallel stratification, and 
the conservation of the total flux (radiative plus convective; the 
convective flux being computed using the local mixing length 
recipe). The radiation field used in the model generation is cal- 
culated by assuming absorption from atoms and molecules by 



opacity sampling at approximately 100000 wavelength points 
over the wavelength range 1300 A-20 /im. 

The code used for calculating the synthetic spectra is 
BSYN v. 7.12 which is based on routines from the MARCS 
code. The atomic line list used in our calculations is com- 
piled from the VALD I database dKupka & Ryabchikova| [l999) 
updated according to Gustafsson et al.l d2008l) . The molecular 
line lists include CO, SiO, TiO, ZrO, VO, O H, H ? Q, CN, C? 
NH, CH, A1 H, SiH, CaH (see refere nces in [G ustafsson et al. 
I2008h. MgH dSkory et al l l2003h . FeH dDulick et al.ll2003l) . and 
CrH dBurrows et al.ll2002l) . Up-to-date dissociation energies and 
partition functions are used. The basic chemical com p ositio n 
adopted is that of the Sun as listed by iGrevesse et a D d2007l) . 
The synthetic spectra were calculated in plane parallel symme- 
try. A constant micro-turbulence velocity of 2 km s _1 is assumed. 
The most important differences to the BT models are the dif- 
ferent opacities sources, the solar abundances, and the fact that 
MARCS is a pure gas phase model with no dust formation. 

Synthetic spectra are calculated for the wavelength region of 
0.50 - 2.53pm, with a resolution of R — 600000. We used a grid 
of MARCS model atmosphere which spans the parameters as 
follows: 

- r eff from 2500 K to 3100K with 100K step (lower tempera- 
tures are not available), 

- logg = 4.5, 5.0, and 5.5 dex, 
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Fig. 4. VLT (NACO) spectra of all the components in J, H and K bands with atomic and molecular lines indicated 
- [M/H] from -0.5 to 0.25 dex with 0.25 dex step. 



3.3. DRIFT 

The DRIFT-PHOENI X mod el atmosphere code has been devel- 
oped b y iDehn et all J2007h ; iHelling et alj d2008bl): IWitte et all 
(120091) . Both the DRIFT code bv lHelling et al.l (l2008al) and the 
BT-Settl code return a consistent dust cloud structure with cor- 
responding opacities and the altitude-dependent depletion and 
redistribution of gas phase abundances, which feed back on both 
the thermodynamical structures and the radiation field. An it- 
eration of these methods allows the determination of stationary 
atmosphere and dust cloud properties and yields the respective 
synthetic spectra. But solely in the case of the DRIFT code, the 
dust formation takes place via the formation of seed particles 
and their subsequent growth or evaporation, solving from top to 
bottom of the atmosphere. The BT-Settl models in contrast solve 
the timescale equations to calculate the depletion of refractory 
elements from bottom to top. 

As in the BT-Settl model, the DRIFT model assumes dirty 
grains mixed according to the composition of each atmospheric 
layer. While the BT-Settl model assumes dirty grains in the 
timescales equations to calculate the growth and settling of the 
grains, it only sums the opacity contributions of each species in 
each layer as for an ensemble of pure grains. The DRIFT model 
instead uses composite optical constants calculated in effective 
mean theory from the composition of the grains dBosch et alJ 



2000). The resulting absorption and scattering properties of the 
grains are therefore different than those of the BT-Settl mod- 
els, possibly producing more opaque clouds. However, since the 
opacities are dominated by atomic and molecular opacities over 
most of the spectral distribution in this spectral type range, the 
impact of those differences are difficult to identify. The largest 
differences between the BT-Dusty, BT-Settl and DRIFT models 
are the differences in the local number density, the size of dust 
grains, as well as their mean composition, which are the direct 
results of the cloud model approach. 

The DRIFT model also considers only seven of the most im- 
portant solids (Ti0 2 , A1 2 3 , Fe, Si0 2 , MgO, MgSi0 3 , Mg 2 Si0 4 ) 
made of six different elements, found to form below an ef- 
fective temperature of 2800 K. Note that the CE, which de- 
termines the composition of the solid species included in 
the BT-Settl model, does not give the same list. Among the 
most important species, forming in the BT-Settl model below 
r eff = 2900 K, are additionally Zr0 2 , CaTi0 3 , CaSi0 3 , Ca 2 Si0 4 , 
Ca 2 Al 2 Si0 7 , Ca 2 MgSi 2 O v , CaMgSi 2 6 , Ti 2 3 , T14O7, Fe, Ni, 
VO, V 2 3 , MgTi0 3 , MgTi 2 5 , MgAl 2 4 , Al 2 Si 2 Oi 3 . The 
DRIFT model includes, similarly to the BT-Settl model, mix- 
ing by convection and overshooting by assuming an exponen- 
tial decrease in mass exc hange frequency in the radiative zone 
(IWoitke & Hellind l2004). However it neglects the contribution 
of the gravity waves included in the BT-Settl model. The model 
code DRIF T-PHOENIX has been discussed in more detail by 
IWitteetal.1 (120091) . 
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We used a grid of DRIFT model atmosphere which spans the 
parameters as follows: 

- from 2200 K to 3 100 K with 100 K step, 

- logg = 4.5,5.0, and5.5dex, 

- [M/H] from -0.5 to 0.5 with 0.5 dex step. 

4. Physical parameters determination and results 

Th e first analysi s of LH S 1070 by spectral synthesis was made 
by iLeinert et all (Il998h using the FOS sp ectra and the pure 
gas phase NextGen model atmospheres dAllard et al.l 1997; 
lHauschildt et all 1 19991) . But the quality of the fits, even for the 
A component was disappointing, especially below 0.72 //m, 
and the T e ff of the B and C components were strongly over- 
estim ated (2700K ) main ly due to the absence of dust forma- 
tion. |Le^e^^2j2000) have therefore used the AMES-Dusty 
models dAllard et al.ll200Th to re-analyse the LHS 1070 system. 
However if the quality of the fits and precision of the for 
the B and C components were clearly improved, this was still 
clearly not the case of the M dwarf primary. The stellar param- 
eters obtained with the assumption of a distance of 8.8 pc are 
summarized: 



- Component A: T eS = 2950 K, logg = 5.3: 

- Component B: T eS = 2400 K, logg = 5.5: 

- Component C: T eB = 2300 K, log g = 5.5: 



:0.2, [M/H] = 0.0 
:0.5, [M/H] = 0.0 
:0.5, [M/H] = 0.0 



In the following, we derive the stellar parameters using more 
recent atmosphere models and spectroscopic informations cov- 
ering both the optical and IR ranges. Metallicity and gravity 
are determined from peculiar spectral features, whereas effective 
temperature and radius are constrained from the overall shape of 
the spectra, following the different steps: (i) a first x 2 minimiza- 
tion is performed on the overall spectra considering effective 
temperature, radius, metallicity, and gravity as free parameters. 
It gives a first guess for the parameter space of each component, 
(ii) we looked for peculiar spectral features that are mainly sen- 
sitive to metallicity (see Sect. 14.2b or gravity (see Sect. 14. 3b to 
refine these two parameters, (iii) we fixed these parameters to 
perform another x 2 minimization and derive effective tempera- 
ture and radius (see Sect. I4.4l i. At each step we checked that the 
resulting value is not sensitive to changes on the value of the 
other parameters. 

Age is estimated from kinematics and rotation. Before en- 
tering the details of our study, we also summarize our results in 
Table [2] Note that thes e values have been obtained by assuming 
a dista nce of 7.72 pc dCosta et alj|2005l) whereas ILeinert et all 
(2000) used the higher value of 8.8 pc. 



4. 1 . Spectral type 

A well-defined spectral classification for the M dwarfs helps in 
the calibration of the temperature of the late type stars and for 
the definition of the end of the main sequence on the HR dia- 
gram. We have derived spectral indices and spectral types for 
all three components. For the early M-dwarf, i.e. the primary, 
we used the classification sch eme based o n the TiO and CaH 
band-strengths, as defined by iReid et all dl995h . For the late 
M-dw a rfs, we also used the PC3 index defined by iLepine et al.l 
(l2003l) : lHawlev et al.l(l2002h : lMartm et aTI(ll999l) . We have com- 
pu ted the spect r al inde x H2O-K in the near IR K-band defined 
by ICovev et al. ( 201 Oh and used the spectral-type vs index re- 
lation from Rojas-Avala et al.l d2010h . The spectral indices and 



corresponding spectral types are given in Table [3] The obtained 
spectral classification for componen t s B an d C is one subclass 
higher as compared to ILeinert et a l. (2000). The spectral types 
obtained in the K-band differs from the optical indices by up to 
two subclasses showing inconsistency on the spectral type ver- 
sus spectral index relations defined from the optical and near-IR 
spectra. Here we adopt the spectral type obtained from optical 
spectral indices. 



Table 3. Spectral index values and derived spectral type com- 
puted from TiO and CaH band-strengths for component A and 
from the PC3 index for components B and C. 



Band 




Indices 


Spectral type 


Component A 


TiO s 




0.211 


M5.5 


CaH 2 




0.281 


M5.5 


CaH, 




0.557 


M6 


H 2 0-K 




0.829 


M6.5 


Adopted value 




M5.5 


Component B 


PC 3 




2.305 


M9.5 


H 2 0-K 




0.791 


M8.5 


Adopted value 




M9.5 


Component C 


PC 3 




2.608 


L0 


H 2 0-K 




0.798 


M8 


Adopted value 




L0 



4.2. Metallicity 

In order to estimate the metallicity of the system, we looked for 
special features in the spectra that are mainly sensitive to the 
metallicity. The main indicator of metallicity is the VO absorp- 
tion band at 7300 A - 7600 A. It is well reproduced at solar 
metallicity for all components, and shown in Fig. [5] for the pri- 
mary with the BT-Settl model. We checked that the same metal- 
licity is found when changing the other parameters. This solar 
metallicity can also be inferred from the Nal, Cal, and H2O fea- 
tures in t he K-band of the primary using the calibration deter- 
mined by Roi as-Avala et al.l (120101) . In the following, we adopt 
this [M/H] = dex value. 



4.3. Gravity 

The surface gravity can be estimated by analyzing the width of 
atomic lines such as the K I and Na I D doublets, as well as the 
relative strength of metal hydrides bands such as those of CaH. 
The K I doublet at 7665 A and 7699 A is a particular useful grav- 
ity discriminant for M stars. Figure|6](left panel) shows the grav- 
ity effects on the strength and pressure broadening of the KI 
lines as modeled by the PHOENIX BT-atmospheres. The over- 
all line strength (central depth and equivalent width) increases 
with gravity as the decreasing ionisation ratio due to the higher 
electron pressure leaves more neutral potassium in the deeper at- 
mosphere. The width of the damping wings in addition increases 
due to the stronger pressure broadening, mainly by H2, He and 
H I collisions. 

Figure [6] (right panel) shows the determination of gravity 
for the component A from the KI doublet using the BT-Settl, 
MARCS, and DRIFT models. The determined log g values are 
given in Table [2] The best agreement is obtained with logg = 
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Table 2. Derived parameters for the LHS 1070 system. The luminosity L is computed from the radius and the effective temperature. 



Component A 



Model 


Tss 


logg 


Radius 


log(L) 




(K) 


(cm s" 2 ) 


CRq) 


(A>) 


BT-Dusty 


2900+100 


5.0 


0.134+0.005 


-2.93+0.090 


BT-Settl 


2900+100 


5.0 


0.134±0.005 


-2.93±0.090 


MARCS 


2900±100 


5.0 


0.136±0.005 


-2.92±0.090 


DRIFT 


2900±100 


5.0 


0.136±0.005 


-2.92±0.031 


isochrones 


2953 


5.2 


0.140 


-2.88 


Component B 


Model 


Tss 


logg 


Radius 


log(L) 




(K) 


(cm s" 2 ) 


C*a) 


(A,) 


BT-Dusty 


2500+100 


5.0 


0.102±0.004 


-3.43+0.105 


BT-Settl 


2500±100 


5.0 


0.102±0.004 


-3.43±0.105 


MARCS 


2600±100 


5.0 


0.098±0.002 


-3.39±0.086 


DRIFT 


2400±100 


5.0 


0.106±0.005 


-3.46±0.044 


isochrones 


2432 


5.3 


0.104 


-3.48 


Component C 


Model 




logg 


Radius 


log(L) 




(K) 


(cm s" 2 ) 


(«q) 


(A.) 


BT-Dusty 


2400+100 


5.0 


0.098±0.002 


-3.53+0.090 


BT-Settl 


2400±100 


5.0 


0.098±0.002 


-3.53+0.090 


MARCS 


2500° ±100 


5.0 


0.100±0.002 


-3.44±0.090 


DRIFT 


2300±100 


5.0 


0.102±0.005 


-3.57±0.034 


isochrones 


2234 


5.3 


0.098 


-3.68 



" Note that this value might be meaningless as there is no cooler available model. 
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Fig. 5. VO band observed in the primary (black) compared to the 
BT-Settl model at 2900 K, log g = 5 .0, R = 0. 1 34fl for different 
metallicities. 

5.0 dex for all the atmospheric models. This value is confirmed 
by the CaH molecular bands (Figure |7). We checked that the 
same metallicity in found when changing the other parameters 
and adopt the value logg = 5.0 dex in the following. For com- 
parison, the gravities inferred from the m asses found from the 
orbits of the system (Seifahr t et alJl2008l) are also listed, which 
are also in good agreement with the values found by fitting. 

4.4. Effective temperature and radius 

We performed a second x 2 minimization by adopting the metal- 
licity and gravity derived in Sect. 14.21 and 14.31 and refined the 
effective temperature and the radius by comparing the overall 
shape of the observed and synthetic spectra. As opposed to the 
studies mentioned previously in which the best fit was found by 
trial and error, in this paper we derive the effective temperature 



and radius by performing a^ 2 minimization technique. For this 
purpose, our approach was to first convolve the synthetic spec- 
trum with a Gaussian kernel at the observed resolution and then 
rebin the outcome with the observation. For each of the observed 
spectra we have calculated the reduced x 2 value by comparing 
these spectra, taking into account their uncertainties estimated 
from the reduction procedures (see Sect. I2.2l i. with the grids of 
synthetic spectra in the wavelength range between 4500 A to 
2.4 yum. We have excluded the spectral region below 4500 A due 
to the low S IN ratio of the observed spectra. The number of data 
points used for the x 2 computation is thus 1487 in the optical 
and 204 in the infrared. 

In a second step, a reduced x 2 ma P nas been obtained for 
each component in the optical and in the infrared as a function 
of temperature and radius. Such a map is shown in Fig. [8] for 
the primary using the BT-Dusty model. The x 2 minimum value 
is given on the left part of the colour bar. The parameter space 
which gives an acceptable solution around the minimum^ 2 val- 
ley is within the white contour, defined by visual inspection. The 
X 2 value at this significance level is indicated along the white 
contour. Error bars are derived from this contour. 

We have identified on the contour maps the possible solu- 
tions in the optical and IR. The adopted values of effective tem- 
perature and radius are the common intersection between the so- 
lutions found in the optical and IR. 

The solutions were finally inspected by comparing it with the 
observed spectra. The same procedure has been used with BT- 
Settl, MARCS and DRIFT model grids. Differences of 100 K to 
200 K in the r e ff determination (see Table |2]i are found for the 
B and C components depending on the model used, whereas all 
models agree on the effective temperature of the primary. 

4.5. Age and mass 

With its high velocity compo nent of W — 43 km/s pe rpendicular 
to the galactic plane based on lBasri & Marcvl d 19951) . LHS 1070 
has been considered as part of the old disk population with an 
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Fig. 8. Left: x 2 ma P computed for component A using the optical spectrum and the BT-Dusty model. Right: same using the IR 
spectrum. All values inside the white contour give an acceptable fit (checked by eye). The adopted common solution to both x 2 
maps is: T e g = 2900 K and a radius of 0.134 R Q . 



age of several Gvrs. iReiners et all d2007bl) refined this estimate 
to about 1 Gyr based on measured rotation velocities of its com- 
ponen ts using a modified Skumanich braking law (Sku manTchl 
1 1 972b . They do not exclude however that the braking law may 
also have to be changed in its absolute time scale, which could 
increase this estimate. 



By an orbital fit, iLeinert et alj d2001l) have computed the 
combined mass of the components B and C and showed that their 
mass is very close to the hydrogen burning minimum mass, in 
good agreement with the masses of 0.0 8 to 0.083MH and 0.079 
to 0.080M o derived by ILeinert et ail (l2000h f r om theoretical 
mass- l uminosity relations {Baraffe et al. 1998; Chabrie r et al.l 
2000J). ISeifahrt et al.l <|2008) constrained the combined mass of 
B an d C to Mb + M c = 0.157 + O.OO9M which is higher 



than I Leinert et al. (2001) because of the refined distance by 
ICosta et al.l ( 2005 ). Recently, an improved fit for the orbit of 
LHS 1070 B and C around each other, and an estimate for the or- 
bit of B and C around A have been performed. The masses of the 
three components are found to be Ma — 0.13 to 0. 1 6M H , Mb = 
0.077 ± 0.005M o , and M c = 0.071 + 0.004M G (IKohler et all 
l20l2h . Here, the values for the primary are uncertain, because the 
wide orbit of this triple system has could not yet been d etermined 
with sufficient accuracy. Finally, Seifa hrt et alJ (120081) have not 
measured individual masses for the B and C components, bu t 
using the mass ratio M C /M B = 0.923 from (iKohler et al.ll2012l) . 
their masses for A, B, C become 0.1 15 + 0.01, 0.082 + 0.01 and 
0.075 + 0.01M Q . 



The interpolations of the NextGen dBaraffe et alJ Il998l) 
isochrone for the primary, and of the AMES -Dusty 
dChabrier et al.l l2000h isochrone for the B and C compo- 
nents, assuming an age of 1 Gyr for the dynamical masses of 
Seifahrt et al. (2008) are shown in Table [2] For masses above 
the hydrogen burning limit these values do not change much 
for larger ages, since in the stellar regime only show small 
evolution effects are seen after the age of 1 Gyr. No evolution 
models have yet been calculated using the BT-Settl models. 
However, only negligible changes are expected with the revision 
of the evolution calculations. Such revision of the interior 
and evolution models is currently being prepared (I. Baraffe, 
Exceter, private communication). 



4.6. Results 

Figures [9] and [10] show the best fit model superimposed to the 
observed optical (FOS) and IR (NICMOS) spectra for all the 
three components using the BT-Settl, BT-Dusty, MARCS, and 
DRIFT models. Note that the observed NICMOS spectra of all 
three components plunge down away from the model predictions 
below 0.85 //m (see Fig. ITOt. whereas the observed FOS spectra 
are correctly represented by the models in the same wavelength 
range (see Fig. [9]). This deviation is due to difficulties with the 
NICMOS data at the very edge of the wavelength range. 

The FOS spectral distribution is better reproduced by the 
models for the primary than for the cooler components B and 
C. The revised opacities and oxygen abundance (among other 
elements) used in the current BT-Settl models allow a signifi- 
cant improvement compared to th e AMES-Dusty and Nex tGen 
models used in previous analysis dLeinert et al.ll9 98. 2000). The 
slope of the spectra is now reproduced over the complete FOS 
spectral distribution, and the strength of molecular bands is re- 
produced in average quite well. However some problems re- 
main which are probably due to uncertain and missing opacity 
sources. Hence, the MgH feature at around 5200 A is too strong 
in all the models while the CaOH band at 5500 A is missing in 
all models. The Nal doublet at around 5900 A is far too strong in 
the models as well as the CaH band at 7000 A. Largest discrep- 
ancies are found around 6000-6400 A for all models. The TiO 
bands around 7055 A as well as the CaH band around 6900 A 
are too strong in all the models. The VO band around 7334 A 
is also visible and is quite well matched by the BT-Settl model. 
The BT-Settl models also differ from the DRIFT models by the 
strength of atomic lines which are deeper in the DRIFT models. 

In the IR NICMOS range, the BT-Settl models fit slightly 
better the primary than the DRIFT models. The MARCS model 
overestimates the flux over most of the spectral range above 1.3 
jum while the DRIFT model shows a slightly different shape of 
the H band peak which is shared by the BT-Settl models in the 
case of the B and C components. But all the models appear over 
luminous in the J and H bands in the case of the B and C com- 
ponents for the selected radius and effective temperature. This is 
also apparent in Figure QT| which shows the comparison of high 
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Fig. 6. Right: BT-Settl and BT-Dusty models for an effective 
temperature of 2900 K and varying log g. The effect of gravity 
and the pressure broadening of the KI doublet is clearly visi- 
ble. The details of the dust treatment only cause negligible dif- 
ferences at this r e ff. Left: KI doublet as observed for the pri- 
mary (black) compared to the BT-Dusty (red), BT-Settl (green), 
DRIFT (blue), and MARCS (brown) models with logg = 5.0. 

resolution IR spectra (NACO) with the best fit of all the four 
models. The change in the Nal and Cal strength as the tempera- 
ture decreases (K-band) is quite well reproduced by the BT-Settl, 
BT-Dusty and DRIFT models. 

As an additional check on the effective temperature determi- 
nation, we also compared the spectrum of the unresolved sys- 
tem in the 3-14 yum range with the BT-Dusty, BT-Settl and 
DRIFT models (MARCS models are not available in this wave- 
length range). The synthetic spectra of the unresolved system are 
computed by adding the individual best fit synthetic spectra. The 
comparison is shown in Fig. [12] where the ISOPHOT spectrum 
is in black, the BT-Dusty model in red, the BT-Settl model in 
green, and the DRIFT model in blue. Star symbols indicates the 
photometry obtained in the IR bands. The overall agreement is 
good except for the observed spectrum above 8.7/mi, where it 
gets quite noisy. 

Relatively small differences distinguish the MARCS, DRIFT 
and BT models in the IR spectral range for the primary. In the 
case of the B and C components, the MARCS models suffer 
clearly of the lack of dust grain scattering which tends to fiat- 
ten out or veil the spectral features in this spectra range. This 
effect is observed in the DRIFT, BT-Settl and BT-Dusty models 
which include dust formation. Differences in the cloud model ap- 



Fig. 7. CaH molecular bands in the spectrum of the primary 
(black) compared to the BT-Settl model at 2900 K, logg = 5.0 
(red) and logg = 5.5 (green). 

proaches explain the differences between the DRIFT, BT-Dusty 
and BT-Settl models. The DRIFT models appear as dusty as the 
2001 AMES-Dusty models with similar effective temperatures 
and su rface gravities for B and C than derived by iLeinert et alj 
(2000). The BT models tends to attribute slightly hotter effec- 
tive temperatures (+100 K) and lower gravities (-0.5 dex) to 
these objects, while the dust-free MARCS models would at- 
tribute them the highest values (+200 K). But judging from the 
overall fits obtained to the NICMOS spectra it appears that nei- 
ther of the models are yet dusty enough to explain the IR spectral 
distribution of B and C. Indeed the over-luminosity of the mod- 
els in the J band could be attributed to missing or to weak veiling 
by dust scattering. 

The Kohler et al. (120 12h mass estimate for the primary re- 
quires a 150 to 300 K higher effective te mperature and up to 
30% larger radius than using the revised Seifa hrt et alJ (2008) 
mass, resulting in an over-prediction of the luminosity of up 
to one magn i tude. In the case of the B and C components, the 
Kohler et al. (2012) mass estimates correspond to effective tem- 
peratures which are 100 to 200 K cooler than obtained in this pa- 
per using the DRIFT models, and in correspondingly larger dis- 
crepancies with the BT-Settl fits. On the other hand, the results 
obtained in this paper for all three components are consistent 
with t he isochrone interpolation for the revised Seifahr t et alj 
(2008]) masses of the A, B and C components of LHS 1070 (see 
Tabled, and there is no evidence that the components may have 
been influenced by their binary nature. 

5. Conclusions 

This paper presents the results from spectral synthesis analysis 
for the LHS 1070 triple system. This system has been exten- 
sively observed from th e optical to the IR, and dynamical masse s 
have been determined (ILeinert et alJl200lHSeifahrt et al.l f2008). 
Therefore, it constitutes a testbed of model atmospheres of low- 
mass stars. Band strength indices are used to measure TiO, CaH 
and PC3 features to classify their spectral type. The components 
are classified as M5.5, M9.5, and L0 dwarfs, and their atmo- 
spheres lie in a temperature range where dust starts to form. We 
have determined the physical parameters T e s, logg, metallicity, 
and radius for the three components of the LHS 1070 system 
by comparing the observed spectra with the synthetic spectra 
computed with the most recent atmospheric models: BT-Dusty, 
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Fig. 9. Optical spectra of all three components. Comparison with model predictions. Black: observed FOS spectra. Green: best fit 
BT-Settl model. Blue: best fit DRIFT model. Brown: best fit MARCS model. The parameters that give the best fit are given in 
Table 2. 
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Fig. 12. Black: ISOPHOT thermal infrared spectra of the un- 
resolved system, with photometric measurements overplotted 
(stars). Red: best fit BT-Dusty model. Green: best fit BT-Settl 
model. Blue: best fit DRIFT model 



BT-Settl, MARCS, and DRIFT. All the models agree for a so- 
lar metallicity for the system. The derived gravity is 5.0 dex and 
agrees within the uncertainties with t he values derived from the 
dynamical mass (ISeifahrtetalJ l2008). 

We found the same value for T e ff for the primary from all 
models while differences of 100K and 200 K are found for 
components B and C depending on the dust density content of 
the model atmosphere used. The revised oxygen abundance by 



lAsplund et ail (120091) and ICaffau et all d2011l) yield significant 
improvements of the BT-Settl fits to the primary compared to 
earlier studies based on the larger valu es of the solar oxygen 
abundance of Grevesse & Sauval (1998). These improvements 
are descri bed inl Allard & Homeierl d20 1 2l) . The even lower abun- 
dances of iGrevesse et alJ (12007 ) used in the MARCS models 
lead to an excess of near-IR flux due to weaker water vapor ab- 
sorption. 

The DRIFT and BT-Settl models differ mainly in their nu- 
merical approach in solving the equations for grain growth, sed- 
imentation and opacities: the DRIFT model solves them from 
the top to the bottom of the atmosphere, while the BT-Settl 
model solves them from the bottom to the top of the atmo- 
sphere. This causes the BT-Settl model to tend to have a deficit of 
grains in the upper atmospheric layers compared to the DRIFT 
model despite an adequate account in both models of supersatu- 
ration effects. Despite these fundamental differences, the result- 
ing grain sizes obtained by the models are quite similar. The 
over-luminosity shared by the models in the J bandpass could 
be indicative of grains of larger sizes and/or more numerous in 
the LHS 1070B and C component atmospheres. The results con- 
firm the Allar d & Homeierl d2012l) findings based on r e ff-color 
constraints. 

The r e ff values found with the DRIFT atmo s phere s for com- 
ponents B and C agree with the iLeinert et al.l (2000) findings. 
Except for the MARCS models which do not include dust treat- 
ment, the models are able to reproduce the observations and 
describe the main features of the visible to IR spectra for all 
three components. This raises the confidence level in the dust- 
modelling approach. 
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Fig. 10. Same as Fig. with NICMOS spectra in the IR. 



However, the calculation of opacities for composite grains 
relies on relatively simple approximations and also does not ac- 
count for possible distributions of grain shapes and structures 
such as porosity. Both models rely on results of radiation hydro- 
dynamical simulations that provide the mixing and overshooting 
which compensates sedimentation effects. One problem is cer- 
tainly that the translation of the resulting radial velocity field 
into a diffusion coefficient is currently uncertain. It is also possi- 
ble that the mixing effects are being currently underestimated by 
local 2D simulations, and that additional mixing is provided by 
other phenomena on larger scales such as global rotation effects. 
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Fig. 11. Black: J (upper panel), H (middle panel), and K (lower panel) NACO spectra of all three components. Comparison with 
model predictions. Green: best fit BT-Settl model. Red: best fit BT-Dusty model. Blue: best fit DRIFT model. Brown: best fit 
MARCS model. The comparison of the BT-Dusty and BT-Settl models together with the observations are shown on the right panels. 
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